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Abstract — In this paper, we have investigated the application of 
non-binary LDPC codes to spatial multiplexing MIMO systems 
with a large number of low power antennas. We demonstrate that 
such large MIMO systems incorporating with low-complexity 
MMSE detector and non-binary LDPC codes can achieve low 
probability of bit error at near MIMO capacity. The new 
proposed non-binary LDPC coded system also performs better 
than other coded large MIMO systems known in the present 
literature. For instance, non-binary LDPC coded BPSK-MIMO 
system with 600 transmit/receive antennas performs within 3.4 
dB from the capacity while the best known turbo coded system 
operates about 9.4 dB away from the capacity. Based on the 
simulation results provided in this paper, the proposed non- 
binary LDPC coded large MIMO system is capable of supporting 
ultra high spectral efficiency at low bit error rate. 

Index Terms — non-binary LDPC code, large MIMO system, 
coded MIMO system, spatial multiplexing, MMSE detection 

L Introduction 

In order to support very high data rate on wireless commu- 
nication channels, it is known that the physical limitation is 
bandwidth since, at the present moment, it is very scarce and 
extremely expensive. One of the possible solutions to meet 
this constraint is the multiple-input multiple-output (MIMO) 
system |1|, |2|. The MIMO system is a transmission system 
that uses multiple antennas at both sides of the communication 
ends. The MIMO system which utilizes spatial multiplexing 
technique (known as BLAST system) fSl, f4l has attracted 
a great deal of attention over the past two decades since it 
provides the significant increase in spectral efficiency without 
additional bandwidth and transmit power in a rich scatter- 
ing environment. The transmitter spatially multiplexes data 
streams and then simultaneously transmits these multiplexed 
data via different transmit antennas. At the receiver side, 
some specific techniques will be employed to demultiplex the 
received data. 

In this paper, we focus our attention to the spatial multi- 
plexing MIMO system that employs tens to hundreds trans- 
mit/receive antennas and refer to such a system as large MIMO 
system. The prominent advantage of large MIMO system is the 
improving in capacity which is proportional to the minimum 
of the number of transmit and receive antennas |2 1. As pointed 
out in the literature O, O, placing a large number of antennas 
is more amenable because the transmitted RF energy can be 
more sharply focused in space. Therefore, we believe that a 
large MIMO system would be practical in the near future to 
support the increasing demand for high data rates in wireless 
communications. 



Concatenating the MIMO system with channel codes is a 
methodology to increase reliability and performance of the 
MIMO system. There has been a tremendous effort to develop 
the coded MIMO systems, which include code-design, the 
invention of soft-output detectors, joint detection and decoding 
techniques and so forth, with the ultimate goal to approach 
MIMO capacity (see e.g. Q, IH, (H). We note that almost 
works have studied the coded MIMO system with less num- 
ber of transmit/receive antennas (e.g. 2 to 8 transmit/receive 
antennas). Therefore, lacking in the literature is a performance 
study of coded large MIMO system. 

To the best of our knowledge, only large MIMO systems 
concatenated with turbo code have been studied in ||3, |[6l, 
1 10|. Although the turbo coded large MIMO system mentioned 
above is very attractive in term of computational complexity 
but there still exists a significant performance gap to the 
MIMO capacity. For example, the gap of turbo coded MIMO 
system with 200 transmit/receive antennas from MIMO ca- 
pacity is more than 7 dB. In this paper, we aim to reduce the 
remaining gap by considering non-binary low-density parity- 
check (NBLDPC) codes El. 

Thank to the superior performance for short and moderate 
codeword lengths, NBLDPC codes have recently received 
an upsurge of research interest from wireless community 
including the application to MIMO channels (121, CS, lfT4l . 
|15|. For the MIMO system with 2 transmit/receive antennas, 
it has been reported in |14| that a regular NBLDPC code 
defined over GF(2^) outperforms both optimized irregular 
binary LDPC code and binary LDPC code defined in the latest 
IEEE standard (the binary LDPC code is specified by parity- 
check matrix defined over GF(2)). However, the NBLDPC 
coded system in (T4\ deploys the optimal maximum a poste- 
rior (MAP) detector to initialize the soft-input for NBLDPC 
decoder. This is impossible for the large MIMO system since 
the complexity of MAP detection grows exponentially with 
the number of transmit antennas and the size of modulation 
constellation. We also note that the performance of NBLDPC 
codes for MIMO systems with suboptimal detector has not 
been reported. 

In this study, we have investigated the application of 
NBLDPC codes to large MIMO systems which utilize the low- 
complexity minimum mean square error (MMSE) detection 
as MIMO detector. Our contributions in this paper can be 
summarized as follows : 

1) The non-binary LDPC coded large MIMO system with 
MMSE detector is proposed in this paper. We also present how 



the soft-output MMSE detector can straightforwardly work 
with NBLDPC decoder. 

2) We provide the simulation results which can be used 
as the benchmark for coded large MIMO systems. We have 
shown that the non-binary LDPC coded system performs best 
among other coded schemes at SNR close to MIMO capacity, 
when a large number of antennas are employed. 

The rest of this paper is organized as follows. We first 
introduce the NBLDPC codes in Section IL In Section III, the 
system model used for all simulations is described. In Section 
IV, we explain the soft-output MMSE detector for NBLDPC 
decoder. In Section V, we present the decoding performance 
of NBLDPC coded large MIMO systems. This paper is closed 
with conclusions. 

II. NBLDPC Codes 

An NBLDPC code C over Galois field GF(2^) is defined 
by the null- space of a sparse P x N parity-check matrix 
A = {aij} defined over GF(2^), for i = 1, . . . ,P and 
j = l,...,7V 

C = {x G GF(2^)^ I Ax"^ = G GF(2^)^}, 

where m > 1 and x = (xi, . . . , xn) is a codeword. The i-th 
parity-check equation can be written as 

anxi + ai2X2 H h aiNXN = 0, 

where a^i , . . . , aiN G GF(2^) are the entries of i-\h row of A. 
The parameter N is the codeword length in symbol. Assuming 
that A is of full rank, the number of information symbols is 
K = N - P diwd the code rate is i? = K/N. 

We note that a non-binary symbol which belongs to GF(2^) 
can be represented by the binary sequence of length m bits. 
For each m, we fix a GF(2^) with a primitive element a 
and its primitive polynomial tt. Once a primitive element a of 
GF(2^) is fixed, each non-binary symbol is given by an m- 
bits representation |16, p. 110]. For example, with a primitive 
element a e GF(2^) such that 7r(a) = a^ + a + 1 = 0, 
each symbol is represented as = (0,0,0), 1 = (1,0,0), 
a = (0,1,0), a'^ = (0,0,1), a^ = (1,1,0), a^ = (0,1,1), 
a^ = (1,1,1) and a^ = (1,0,1). Let L{x) be the binary 
representation of x G GF(2^). For the above example, we 
can write L{x = a^) = (1,1,0). Thus, each coded symbol 
Xi e GF(2^),Vz G {1,...,A^} of a non-binary codeword 
represents m bits. We also denote n = mN and k = mK as 
the codeword length and information length in bit, respectively. 

An NBLDPC code is (dv, <^c) -regular if the parity-check 
matrix of the code has constant column weight d^ and row 
weight dc. The parity-check matrix A can be represented by 
a Tanner graph with variable and check nodes (TTl p. 75]. 
The belief propagation (BP) algorithm for NBLDPC decoder 
ifTTIl exchanges the probability vector of length 2^ between 
variable nodes and check nodes of the Tanner graph at each 
iteration round I. 

In this paper, only is {d^ = 2, dc) -regular NBLDPC code 
defined over GF(2^) considered due to the following reasons 

1) The process to optimize parity-check matrix A is not 
required since it is empirically known as the best performing 



code especially for short code length. Moreover, the NBLDPC 
code with d^ = 2 can be encoded in linear time 1 18]. 

2) The high decoding complexity of NBLDPC decoder can 
be compensated since the A is very sparse 1 191 . 

3) As stated in the introduction, NBLDPC code defined 
over GF(2^) seems to offer the best performance for MIMO 
system with MAP detection |[T4l . We therefore expect the 
excellent performance of NBLDPC codes when applying to 
large MIMO systems. 

4) It is empirically shown that the application of this code 
to higher order modulation is outstanding |20, p. 32]. 

It is not overly exaggerated to state that we intend to apply 
the simple, low-complexity, high-performance channel code to 
large MIMO systems. 

III. System Model 

We adopt the conventional notation to denote the MIMO 
system with Nt transmit antennas and Nr receive antennas as 
Nt X Nr MIMO system. Let A^ be the complex modulation 
constellation of size M = 2^ where p represents bit(s) 
per modulated symbol. In this study, each antenna uses the 
same modulation scheme and the mapping is a Gray-labelled 
constellation. 

Figure [T] shows the spatial multiplexing MIMO system 
concatenated with an NBLDPC code of rate R defined over 
GF(28). 
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Fig. 1. NBLDPC coded MIMO system. Every 8 bits is mapped to a symbol 
in GF(28). 



At the transmitter side, a bit to symbol mapper maps a group 
of 8 information bits to a non-binary symbol in GF(2^). The 
stream of K non-binary symbols is encoded into a codeword 
of length N symbols through an NBLDPC encoder. Each 
coded symbol in GF(2^) is then mapped to a group of q = 
8/p modulated symbols by a constellation mapper. At each 
time instant (each channel use), the transmitter simultaneously 
sends Nt = Ktq modulated symbols in parallel through Nt 
transmit antennas where Kt is a number of coded symbols per 
each transmission. Let s = [si, S2, . . . , ^atJ^ G C^* be the 
transmit signal vector. Each entry s^, Vz G {1, . . . , A/"^} taken 
from A^ must satisfy the component-wise energy constraint 
^[ll^ilP] = Es/Nt where Eg is the total transmitted power 
and E[-] denotes the expectation. With this energy constraint, 
a large number of transmit antennas imply low power con- 
sumption per each transmit antenna. 

Consider a 16 x 16 MIMO system with QPSK modulation 
(M = 4 and p = 2) as an example. After encoding, each 



coded symbol in GF(2^) is mapped to g = 4 modulated 
symbols. At each time instant, the transmitter collects 16 
modulated symbols mapping from Kf = 4 coded symbols. 
These 16 modulated symbols are multiplexed and sent through 
16 transmit antennas. 

Considering a baseband discrete time model for a flat fading 
MIMO channel, the received vector y = [?/i, ?/2, • • • , UNr-V ^ 
C^^ of the spatial multiplexing Nt x Nr MIMO system is 
given by [5] 

y = Hs + n. (1) 

The matrix H = [H1H2 . . . Hat,] G C^-x^* denotes the 
channel fading matrix whose entry hkj is assumed to be 
complex Gaussian random variable with zero mean and unit 
variance E[||/i/cj|p] = 1. The vector n = [ni, n2, . . . , n^^]^ G 
(^Nr- jg ^ noise vector whose entry is a complex white Gaussian 
noise with zero mean and variance a^ per real component. 
The MIMO detector performs detection and produces the 
prior probabilities (soft-output) for NBLDPC decoder. After 
all N variable nodes are initialized, the NBLDPC decoder 
performs decoding process and provides the estimated non- 
bianry symbols (hard output). These estimated symbols are 
finally demapped to a sequence of estimated information bits. 
In this paper, the channel matrix H is assumed to be known at 
the receiver and we only focus on the square channel matrix, 
i.e. Nt = Nr. 

Since each entry of H has unit variance, the average signal 
energy per receive antenna is Es. We follow the convention 
that No/2 = a'^ to define the signal to noise ratio. In this 
setting, the average signal to noise ratio (SNR) per receive 
antenna, denoted by 7, is given by fl} 

7=# = #^. (2) 



A^n 



2ar 



The spectral efficiency (transmitted information rate) of 
coded MIMO system with spatial multiplexing technique is 
pRNf |7|. With perfect H at the receiver side, ergodic MIMO 
capacity is given by [2J 



C = E [log2 det (I^, + (j/Nt) HH^)] 



(3) 



where the superscript H denotes the Hermitian transpose 
operator, det denotes the determinant and In^ is the identity 
matrix of size NrXNr. The MIMO capacity defined in ^ can 
be evaluated via the Monte-Carlo simulation. Both of spectral 
efficiency and capacity are measured in bits/sec/Hz (bps/Hz). 

IV. Soft-Output MMSE Detector 
FOR NBLDPC Decoder 

One of the major problems that prohibits large MIMO sys- 
tems from practical implementation is the high computational 
complexity of MIMO detector. Even with modern and high- 
speed circuits, the optimal MAP or ML (maximum likelihood) 
detector is still infeasible. Therefore, many works have re- 
cently focused on inventing the low-complexity detectors to 
enable the use of large scale multiple antennas I2TII . 

In this paper, we consider the MMSE detector which is a 
famous class of linear detection since it has low computational 
complexity. Based on the complexity analysis provided in 



(5], the MMSE detector has an average per-bit complexity of 
0{NtNr) which is lower than the well-known ZF-SIC (zero- 
forcing with successive interference cancellation) detector 
whose average per-bit complexity is 0{N^Nr). Several near- 
optimal MIMO detection algorithms such as sphere detection 
or lattice reduction-aided detection also exhibit high compu- 
tational complexity which is much greater than 0{NtNr) 10. 
Let us introduce the MMSE detector and the soft-output 
generation for NBLDPC decoder. Following |22|, a detection 
estimate Sk of the transmitted symbol on k-th antenna by 
multiplying i/k with the MMSE weight matrix W/e 



Sk =W]^yk, 



where the MMSE weight matrix W^ is of the form 

-1 



W^ = 



Nc 







Es/Nt 



>-N^ 



HH 



Hk' 



(4) 



(5) 



This MMSE weight matrix W/e is chosen so as to minimize 
the mean square error between the transmitted symbol Sk and 
5 /e. It is important to note that direct computation of the matrix 
inverse in ([5| can be avoided by using Sherman-Morrison- 
Woodbury formula EH p. 50] 

The estimation s^ can be approximated as the output of 
an equivalent additive white Gaussian noise (AWGN) channel 
L24J 

Sk= iJikSk^ Zk, (6) 

where (ik = W^H/g and z^ is a zero-mean complex Gaussian 



variable with variance e^ = 



jfi/iik ~ l^\)' Based on this 



approximation, the probability of Sk conditioned on s G A^ 
is as follows 



Pr(4 \s) 



rexp 



1 



Sk - l^kS 



(7) 



where || • |p denotes the squared Euclidean norm and k, is the 
normalized constant such that XI^^^m Pr {sk \ s) = 1. 

Let pv \x) denote the input of belief propagation (BP) 
algorithm which is used as NBLDPC decoder where v = 
1, 2, . . . , A" and the index in superscript represents the iteration 
round of BP algorithm. The pv ^ (x) is the probability that v-th 
coded symbol is likely to be x G GF(2^). We assume that k-th 
transmit antenna to {k -\- q — l)-th transmit antenna are used 
to send the v-th coded symbol (any coded symbol is mapped 
to q modulated symbols). We mention that the modulated or 
constellation symbol s^ G A^ can be demapped to p bits 
according to the label of constellation. For each x G GF(2^), 
we need to collect q modulated symbols to represent qp = S 
bits. The generation of soft-output from MMSE detector for 
v-th variable node of BP algorithm is as follows. 



g-i 



pi'\x) 



Y[ Pr {h+i 



ofc + i 



(8) 



i=0 



where Sk-\-i is the estimation on the {k -\- i)-th receive an- 

^ and the eight bits ordered sequence of 

fe+i . . . ok+q-i^ must be equal to the binary represen- 



tenna, 5^+* G 



{s^,s-'%--'S- 

tation L(x). According to ([8]), 2^{q — 1) real multiplications 

for each coded symbols are needed to calculate pv (x) for 



all X e GF(2^). We note that the computational complexity 
of generating soft-output for NBLDPC decoder is rather low 
when comparing to the computation of MMSE matrix. 

For the sake of completeness, we now introduce the BP 
algorithm for NBLDPC decoder. The BP algorithm mainly 
consists of 4 parts which can be described as follows (in this 
paper we use m = 8). 

Initialization : We set the iteration round i = and define 
the maximum iteration ^max- For each variable node v, pv (x) 
is computed from (l8l. Each variable node sends the initial 
message pic~ = Pv ^ I^^"^ to each adjacent check node c 
where c = 1, . . . ,P. 

Check to Variable : For each check node c, let dc be the set 
of adjacent variable nodes of c. The check node c sends the 
following message p\ 
V edc 



' ^ G M^"" to each adjacent variable node 



Fcv 



P^v}{<^v}x) for a; G GF(2"'), 



P, 



(^+1) 



ix) 



where pi ® P2 ^ 



= pii-^^\acvx) for X e GF(2^), 
i^"^ is the convolution of pi G I 



and 



P2 e 



which can be expressed as 



(Pl ^P2){x) 



^ p,{y)p2{z)forxeGF{2^). 



y,zeGFi2'^) 
x=y-\-z 



The convolution appeared above can be efficiently calculated 
via FFT and IFFT |25| . Increment the iteration round as ^ := 
^ + 1. 

Variable to Check : For each variable node v = 1,2,..., A/", 
let dy be the set of adjacent check nodes of v. The message 



Pv 



G GF(2^) sent from v to c is computed by 



J^)/^A_eJ0). 



.W- 



P^c^(x) = ipr{x) Cc'eaAlcl^c^U^) for X e GF(2-), 
where ^ is the normalized constant so that 

Z^a^GGF(2"^) P'^c [X) — 1. 

Tentative Decision : The tentative decision of the v-th symbol 
is given by 

x^^^ = argmaxp^°)(x) Y[ pi^J{x). 

a:GGF(2-) ^^^^ 

The algorithm stops when the maximum iteration ^^ax is 
reached or Ax = G GF(2^)^. Otherwise repeat the latter 
3 decoding steps. 

V. Simulation Results 

In this section, the performance of NBLDPC coded MIMO 
system with MMSE detector (NBLDPC coded system for 
simplicity) is presented. The maximum iteration of NBLDPC 
decoder ^^ax is set to 100 for all simulation results. 

In Fig. [2j we plot the bit error rate (BER) performance of 
NBLDPC coded 16 x 16 MIMO system with QPSK modula- 
tion (M = 4). The turbo coded STBC (space-time block code) 
MIMO systems with LAS (likelihood ascent search) detector 
161 are chosen for comparison purpose. For the 16 x 16 MIMO 
system with QPSK modulation, these turbo coded systems 
perform within about 4 dB from the capacity which is known 



as the best performance founded in the literature. By using 
R = 1/3 and R = 1/2 NBLDPC codes with n = 3456 bits, 
NBLDPC coded systems outperform turbo coded systems by 
approximately 0.5 dB. We also show the BER performance of 
short length NBLDPC codes n = 864. As we expected, the 
use of shorter length n = 864 codes result in degraded BER 
performance. In term of nearness-to-capacity, the performance 
of NBLDPC coded systems is 3.5 dB away from MIMO 
capacity. We mention that the MIMO transmission can be 
classified into two main categories : 1) spatial multiplexing 
for higher data rate and 2) space-time coding for higher 
transmission quality. We, however, compare the turbo coded 
system (61 which is space-time coding and the proposed 
NBLDPC coded system which is spatial multiplexing since 
both have the same spectral efficiency (also MIMO capacity). 
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Fig. 2. Bit error rate curve of coded 16 x 16 MIMO systems with QPSK 
modulation. The spectral efficiencies of 10.6 and 16 bps/Hz are obtained from 
MIMO system with channel codes of i^ = 1/3 and R = 1/2, respectively. 



The BER performance of NBLDPC coded 200 x 200 MIMO 
systems with BPSK modulation (M = 2) is shown in Fig. Is] 
For 200 X 200 MIMO system, the best performing scheme 
which can be founded in the Hterature is the R = 1/2 turbo 
coded spatial multiplexing MIMO systems with the MMSE- 
LAS detector |5|.The MMSE-LAS detection algorithm uses 
the MMSE detection to initialize the algorithm. Therefore, the 
overall computational complexity of MMSE-LAS detector is 
greater than that of MMSE detector. The BER performance of 
this turbo coded system is away from MIMO capacity by 7.5 
dB. It is clearly seen from the figure that R = 1/2 NBLDPC 
coded system with n = 2400 bits significantly outperforms 
turbo coded system by about 4 dB. Both R = 1/3 and R = 
1/2 NBLDPC coded systems perform within just 3.5 dB from 
the MIMO capacity 

Figure |4]presents the simulation results for larger dimension, 
600 X 600 MIMO system. It can be observed that R = l/3 and 
R = 1/2 NBLDPC coded systems operate within 3.4 and 3.6 
dB, respectively, from the corresponding MIMO capacities. 
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Fig. 3. Bit error rate curve of coded 200 x 200 MIMO systems with BPSK 
modulation. The spectral efficiencies of 66.6 and 100 bps/Hz are obtained from 
MIMO system with channel codes of R = 1/3 and R = 1/2 respectively. 



The turbo coded system with MMSE-LAS detector in |T| is 
again used as the competitor. The performance of this turbo 
coded system over 600 x 600 MIMO system is away from 
the capacity by about 9.4 dB. It is obviously seen that R = 
1/3 turbo coded system underperforms R = 1/3 NBLDPC 
coded system by more than 6 dB. More interestingly, R = 
1/2 NBLDPC coded system absolutely outperforms R= 1/3 
turbo coded system by about 3.4 dB. Therefore, the NBLDPC 
coded system outperforms the turbo coded system both in term 
of performance and spectral efficiency. 
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Fig. 4. Bit error rate curve of coded 600 x 600 MIMO systems with BPSK 
modulation. The spectral efficiencies of 200 and 300 bps/Hz are obtained from 
MIMO system with channel codes ofR= 1/3 and R = 1/2 respectively. 

We also investigate the BER performance of coded large 
MIMO system with higher-order modulation. The coded per- 
formance of 600 X 600 MIMO system with 16-QAM is shown 



in Fig. [5] and [6] From both figures, the following observations 
can be listed as follows : 

• R = l/3 mid R = l/2 NBLDPC coded systems with n = 
2400 bits operate within 6 dB and 8 dB, respectively, from 
MIMO capacity. 

• R = 1/3 and R = 1/2 turbo coded systems |5| operate 
very far (more than 15 dB) from MIMO capacity. It is clearly 
seen that the NBLDPC coded systems indeed beat turbo coded 
systems. 

• By increasing the length of code from n = 2400 bits to 
n = 28800 bits, the coding gain about 0.7 dB can be obtained 
from R=l/3 NBLDPC coded system. 

• For moderate codeword lengths (n = 2400 bits) and the 
same MMSE detector, R = 1/3 NBLDPC coded systems 
outperform both the optimized and regular binary LDPC 
(BLDPC) codes. The coding gains obtained over optimized 
and regular binary LDPC codes are about 0.8 and 2 dB 
respectively. In Fig. [6j the regular binary LDPC code has 
column weight 4 and row weight 6 whereas the optimized 
binary LDPC code is taken from Table. Ill in |26|. The 
maximum degree of variable node of the selected optimized 
binary LDPC code is 16. This optimized LDPC code can 
asymptotically performs very close to the capacity of single 
input single output Rayleigh fading channel (within 0.19 dB). 

• Another advantage of using NBLDPC code which can be 
seen from this figure is the excellent frame error rate (FER). 
Comparing with optimized BLDPC code, sharp decrease in the 
FER curve is obtained. Although the BER curve of optimized 
BLDPC code is good but its corresponding FER curve is quite 
bad. 

For binary LDPC coded MIMO system with MMSE detec- 
tor, the soft output from detector is slightly different from 
NBLDPC coded system. We note that the input of binary 
LDPC decoder is the probability of size two which can be 
represented by log likelihood ratio (LLR). The LLR for Ith 
bit (V/ G {1, . . . ,p}) of the symbol Sk transmitted from kth 
antenna is given as 



In 






where A^^ is the set of modulated symbols whose binary 
labelling of Ith bit is 1 (A^ is similarly interpreted) and 
Pr (sfc I s) is computed from ([t]). It is known that the per- 
formance of BLDPC coded system can be enhanced by 
adding the number of iterations between decoder and detector 
but, for large MIMO systems, this will greatly increase the 
computational complexity. Thus, the simulation for BLDPC 
coded system with joint detection and decoding is not carried 
out. 

Although the MMSE detector has very poor detecting 
performance for MIMO system with Nf = Nr but in this 
section we have demonstrated that the application of this 
detector to NBLDPC coded system is somewhat excellent. 
The uncoded performance of MMSE detector is extremely 
poor comparing to the optimal one. The MMSE detector has 
a diversity order of Nr — Nf -\- 1 while the optimal ML 
detector has a diversity order of Nr |5|. The MMSE detector 
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MIMO detection such as the one described in (TT, p. 358] leads 
to seriously degraded performance (not shown here). 
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Fig. 5. Bit error rate performances of coded 600 x 600 MIMO systems with 
16-QAM modulation. The spectral efficiencies of 800 and 1200 bps/Hz are 
obtained from MIMO system with channel codes of i^ = 1/3 and R = 1/2 
respectively. 
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Fig. 6. Performance comparison between R = 1/3 NBLPDC code and R = 
1/3 BLDPC codes on 600 x 600 MIMO systems with 16-QAM modulation 
at n = 2400 bits. The solid curves represent the BER performance while the 
dashed curves are the corresponding FER performance. 



thus achieve only first order diversity for MIMO system with 
Nt = Nr. Therefore, One may be surprised why MMSE 
detector yields such good coded performance. Let us take a 
closer look at the uncoded results shown in Fig. [7] It can 
be observed that the uncoded performance of MMSE detector 
for 600 X 600 and 16 x 16 MIMO systems are very similar to 
those of near-optimal MIMO detectors for low SNR region. 
Noticeably, those regions are the operating regions for coded 
systems. By trial and error, we suspect that the key component 
to achieve good performance from NBLDPC coded system 
is not a detection part but a soft-output generation. We note 
that other methodologies for generating soft-output from linear 
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Fig. 7. Uncoded performances for 600 x 600 and 16 x 16 MIMO 
systems. The results of MMSE-LAS detection |5|, Tabu search |28|, sphere 
detection |28 1, Belief propagation (BP) based detection |29 1, MMSE-OSIC^ 
detection [30 1 which have very good performance in large MIMO systems are 
included for comparative purpose. The operating regions for coded system are 
highlighted by ellipses. 

Based on the simulations provided in this section, we list 
our major results as follows : 

1) The BER performances obtained by NBLDPC coded 
large MIMO systems are now the best known performance 
and very useful for further investigations. 

2) In term of BER performance, the NBLDPC coded sys- 
tems outperform turbo coded systems and also binary LDPC 
coded systems. The provided results reveal that the turbo code 
may not be a powerful channel code for large MIMO systems 
especially when the higher order modulation is adopted. 

3) The application of MMSE detector to coded large MIMO 
systems is very promising if the operating region is near the 
MIMO capacity. 

VI. Conclusion 

In this paper, the NBLDPC coded large MIMO systems 
are studied. The low-complexity MMSE detector is employed 
to provide the soft-input for NBLDPC decoder. We have 
demonstrated that the proposed NBLDPC coded system can 
significantly decrease the remaining gap from MIMO capacity 
which is previously obtained from the best known turbo coded 
systems. By using moderate length NBLDPC codes (a few 
thousand bits), the proposed coded system can perform near 
MIMO capacity which is closest than ever. We therefore 
conclude that the NBLDPC coded large MIMO system can 
be one of the best choices to achieve both the excellent BER 
performance and the ultra high spectral efficiency. 
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